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We study theoretically the proximity effect of a one-dimensional metallic quantum wire (in the absence of 
spin-orbit interaction) lying on top of an unconventional superconductor Three different material classes are 
considered as a substrate: (i) a chiral superconductor in class D with broken time-reversal symmetry; a class 
Dili superconductor (ii) with and (iii) without a nontrivial Z 2 number Interestingly, we find degenerate zero 
energy Majorana bound states at both ends of the wire for all three cases. They are unstable against spin-orbit 
interaction in case (i) while they are topologically protected by time-reversal symmetry in cases (ii) and (iii). 
Remarkably, we show that non-local spin correlations between the two ends of the wire can be simply controlled 
by a gate potential in our setup. 

PACS numbers: 74.45+c, 74.20.Rp 



Introduction - Proximity effects of superconductors and 
normal metals as well as semiconductors have been a sub- 
ject of continued interest due to the rich physical phenomena 
of these hybrid systems. Especially, the topological nature of 
the superconducting proximity effect is an issue of current re- 
search activities [1-4] stimulated by the possible occurrence 
of Majorana fermions in solid state systems [5-11]. For ex- 
ample, the interface between an s-wave superconductor and 
the surface state of a three-dimensional topological insula- 
tor (TI) is predicted to offer a platform to realize Majorana 
fermions [1]. The underlying idea here is holographic princi- 
ple that entails the reduction of electronic degrees of freedom, 
i.e., electron fractionalization, at the surface of a topologically 
non- trivial bulk state. Furthermore, the proximity effect of 
a TI to unconventional superconductors has also been stud- 
ied theoretically [12]. Another (highly interesting) proposal 
is to use semiconductors with Rashba spin-orbit interaction 
in combination with s-wave superconductors for this purpose. 
In the presence of a sufficiently strong magnetic field, an in- 
verted gap opens at the T-point in this system and topological 
superconductivity as well as Majorana fermions may appear 
when the Fermi energy lies inside the gap [13-16]. A recent 
experiment observed the zero-bias conductance peak at the 
edge of an InSb quantum wire (QW) coupled to an s-wave su- 
perconducting substrate, which might be the first experimen- 
tal observation of a Majorana fermion [17]. It has also been 
shown theoretically that a one-dimensional Rashba quantum 
wire coupled to a d-wave superconductor hosts doubly degen- 
erate Majorana bound states (MBSs) at the edge [18]. Besides 
these distinct works, there are various other research activities 
on the basis of QWs (for example Refs. [16, 19-25]). In all the 
cases listed above, spin-orbit interaction and/or the Zeeman 
energy in the QW is essential. Therefore, an open question is 



if zero energy Majorana bound states are also realizable in a 
one-dimensional system in the absence of these interactions, 
which opens new opportunities as we will show below. 

In this paper, we investigate electronic states caused by 
the proximity effect between a metallic QW and a two- 
dimensional (2D) unconventional superconductor to search 
for zero energy MBSs at the ends of the QW. We discover 
that the resulting MBSs are doubly degenerate, i.e., charac- 
terized by a spin degree of freedom. This leads to nonlocal 
spin-correlations between the two ends of the QW which can 
be manipulated by all-electric means. For a substrate in sym- 
metry class D [26], i.e., in the absence of time reversal sym- 
metry (TRS) the degenerate end states can be gapped out by 
switching on TRS preserving local imperfections that couple 
opposite spin. This is reflected in the fact these systems are 
characterized by a trivial Z2 topological invariant. In contrast, 
for a TRS preserving substrate in symmetry class Dili, we find 
helical MBS pairs that are topologically protected by TRS. 

Model- Our setup is shown schematically in Fig. 1 and the 
model Hamiltonian given by 
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FIG. 1. One dimensional metallic quantum wire (QW) on top of 
a substrate of an unconventional two dimensional SC. The wire is 
located along the x axis. Dashed lines illustrate tunneling between 
the QW and the SC, and solid spheres represent the sites of our tight- 
binding model discussed in the text. 



where i^ w ire describes the QW put along x-direction on the 
surface of the superconductor described by Hsc which in- 
finitely spreads in the xy-plane. These two are connected by a 
tunneling Hamiltonian H t with the spin-independent constant 
matrix element c. We consider the following three cases for 
the order parameter matrix A: (i) A = ia y Aa z (k x + ik y ) 
for the chiral superconductor of class D with the topological 
invariant Z = 1, and (ii) A = ia y A k • a for the time rever- 
sal symmetric class Dili superconductor with the topological 
invariant Z2 = 1. Here, the a denotes Pauli matrices in spin 
space. As the third case (iii), we consider a system where 
two layers of a 2D superconductor of case (ii) are coupled to 
each other by an interlayer transfer integral, which is a topo- 
logical^ trivial case with Z2 = 0. The model (i) is relevant 
to Sr2RuC>4 where the chiral p + ip superconductivity with 
broken time-reversal symmetry is believed to be realized [27] . 
The model (ii) is the 2D analogue of the He B phase, and 
is continuously connected to the helical non-centrosymmetric 
topological superconductor with Rashba spin-orbit interac- 
tion [28-30]. We present the calculational steps and results 
mainly for case (i) below, but it is straightforward to apply 
them to the other two cases (ii) and (iii). Our main purpose 
is to investigate the electronic states in the one-dimensional 
superconductivity formed in this way and to discuss its topo- 
logical nature as well as resulting Majorana bound states. 

Effective Green's function - First, we derive the effective 
Green's function of the QW by integrating over the supercon- 
ductor 



I(ioj n ,k x ) = J dk u 



wire (k x )r z - c I(icj,k x ), (5) 



(6) 



where r denotes Pauli matrices in the Nambu space. As long 
as there is neither spin-orbit interaction in the wire nor spin- 
dependent tunneling between the wire and the superconduc- 
tor, theinduced superconducting order parameter is also spin- 
triplet and can be decoupled into two sectors a x = ±1. The 
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FIG. 2. Spectral function of the electronic excitations. The insets 
show the schematic structure of the spectrum; Blue: continuum of 
the substrate superconductor, Red: ID states of the QW, and Green: 
edge states of the 2D SC induced by the coupling to the QW. In (a) 
the states in the QW are located outside the superconducting gap, 
while in (b) these states are within the superconducting gap and an 
SC order parameter is induced in the QW. (c) is the case where the 
coupling is strong enough that the bonding and anti-bonding states 
act as a potential barrier and edge states in the 2D superconductor 
emerge, (d) shows that no zero energy states are formed in an s-wave 
superconductor under the same conditions as (c). 



integration can be analytically performed and we obtain 

G~ff(iu n , k x ) = iu n (l - 2c 2 h) - h{k x ) • r, 
h(k x ) = (2c 2 Ak x ha x , 0, & + 2c 2 f k c h + 2c 2 J 2 ) , 
m s [A^ 1 arctan («/A+) — Al 1 arctan (ac/A_)] 

1 ~ 27r v / ra 2 A 4 - 2ra s A 2 /i s - uj* 

X- arctan (k/X-) — A + arctan (k,/\+) 

2 ~ 47r v / ra 2 A 4 - 2m s A 2 /i s - J n ' 

is a momentum cut-off and X± 
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where k 

(kl + 2m 2 A 2 - 2m s p s ± 2m sV / m 2 A 4 - 2m s A 2 p s - u\ 

£k — k 2 /2m — p, ^ c likewise. In the following, we look at 
the a x = 1 sector, which is equivalent to a spinless p-wave 
superconductor. From Eq. (8), we can then calculate the 
spectral function A(uj, k x ) of the electron as 

A[u,k x ) = - ^3 [Tr [r z Gfv(oJ,k x )]] , (11) 

where G^(lj, k x ) is the retarded Green's function obtained 
from G e ff (iu) n , k x ) via the analytic continuation iuj n — >> uj + 
iS, where S is a infinitesimally small positive number. In 
Fig. 2, we plot Eq. (11) with m = rn s = 1.0, p s = 1.0, A = 
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0.5, and several values of /i and c. The insets show the 
schematic structure of the spectrum. The red lines represent 
states in the wire without and with an induced superconduct- 
ing gap, in cases (a) and (b) respectively. The blue shaded area 
shows the continuum of the bulk states in the superconductor, 
which arises in the effective Green's function of the QW as the 
self-energy (Eqs. (5) and (6)). The green lines are edge states 
in the 2D SC that will be further discussed below. Figures 2(a) 
and (b) correspond to the cases where the band of the QW is 
outside and inside of the superconducting gap of the substrate 
superconductor, respectively. We can see that in Fig. 2(b) a 
finite gap is induced by the proximity effect. An interesting 
situation arises when the magnitude of the transfer integral 
related to the tunneling c is strong enough that bonding and 
anti-bonding states are formed at the interface of the wire and 
the superconductor. The energy of those states is pushed away 
from the low energy regime, and they act as a potential barrier 
to other states. Therefore, a boundary is effectively induced in 
the substrate superconductor. It is known that edge modes ap- 
pear at the boundary of 2D topological superconductors. The 
energy of the midgap states around k x = approach zero as 
the magnitude of c becomes larger (see Fig. 2(c)). This shows 
the formation of edge channel that are gapped out because the 
edge modes located at both sides of the wire weakly interact 
with each other. In Fig. 2(d) we show the spectral function for 
an 5-wave superconductor substrate for comparison. There 
are no states corresponding to the edge channel in case (c). 

Topological nature - Now, we focus on the topological na- 
ture of the QW system. For this purpose, we use the method 
for the calculation of the topological number in terms of 
Green's function at zero frequency [31, 32]. We derive an 
effective Hamiltonian from Eq. (5) on the basis of H e ^{k x ) = 
-G~fi(iuj n = 0,k x ) = h(k x ) ■ r for a given cr x = ±1. 
This 2x2 Hamiltonian is similar to that of Kitaev's origi- 
nal proposal of a ID topological superconductor [33]. It 
is known that the system is in the topological phase when 
h(k x ) = h/\h\ connects antipodal points of the unit sphere 
as k x is varied from the center to the boundary of the Bril- 
louin zone [4]. Therefore, it is easy to verify that when 
fi > 2c 2 (-fiJ 1 (k x = 0) + I 2 (k x = 0)) = /i*, the QW be- 
comes two copies of Kitaev's ID topological superconductor 
with Z 2 = 1. The topological nature of the ID superconductor 
manifests itself in bound states at the ends of the QW, which 
will benumerically studied in the following. However, in case 
(i), the two identical copies of the Z 2 -nontrivial Kitaev model 
form a trivial composite system. Physically, this implies that 
the degenerate end states can be gapped out by switching on 
terms that couple the opposite spin sectors as we confirmed 
by numerical calculations. 

Numerical study of a tight-binding model and Majorana 
bound states.- We study a tight-binding model by numeri- 
cal calculations. We change our continuum model to a lattice 
model, i.e., k —> sin k and k 2 —> 2(1 — cos k) in Eq. (8), and 
construct the corresponding tight-binding Hamiltonian. Then, 
we calculate the energy spectrum in the system with open 
boundary conditions as depicted in Fig. 1. Figure 3(a) shows 



the topological phase transition as a function of the chemi- 
cal potential in the QW. In the strong coupling regime, where 
H < /i*, there are no states in the induced superconducting 
gap, while in the weak coupling regime, where /i > p* there 
are zero energy states. It can be checked that four-fold degen- 
erate zero energy states are formed, two at each end. Fig. 3(b) 
shows the probability distributions of one of them. Red solid 
circles represent the weight in the QW, and we can see the 
state is localized at the ends of the wire. In sharp contrast to 
the Kitaev model, we have two spin sectors in the QW which 
are degenerate. This fact attaches a spin degree of freedom to 
the MBSs. However, in the presence of spin-orbit interaction, 
which mixes the a x = ±1 sectors, the degeneracy is lifted and 
the Majorana fermions will be pushed away from zero energy. 

In the following, we briefly discuss the other two possibili- 
ties (ii) and (iii) for the SC substrate. We have confirmed that 
a similar analysis to case (i) outlined above applies to these 
cases, and we have found that the doubly degenerate zero en- 
ergy Majorana bound states appear when the chemical poten- 
tial p of the QW is positive and large enough for the topolog- 
ical phase transition to occur. In cases (ii) and (iii), the spin 
degenerate pair of MBS stems from a non-trivial Z 2 invariant 
for the composite system in symmetry class Dili. Therefore, 
these end states are topologically protected by TRS and cannot 
be gapped out by switching on spin-orbit interaction as long 
as the bulk gap is maintained. This is an essential difference 
to case (i) where the conservation of cr x = ±1 was needed 
in addition to the generic symmetries of the model to obtain 
zero energy MBSs, while TRS protects the doubly degenerate 
MBSs in cases (ii) and (iii). This is particularly remarkable for 
case (iii) where the substrate 2D superconductor is originally 
in a topologically trivial phase. 

Local and non-local fermions - There have been many the- 
oretical proposals to use chiral MBSs for topological quan- 
tum computing [34-36]. Chiral MBSs are difficult to de- 
tect and to manipulate because they only couple in a lim- 
ited way to their environment. This is different in our pro- 
posal. On the one hand, the doubly degenerate MBSs dis- 
cussed above are more susceptible regarding their coupling to 
the environment. On the other hand, they can show novel re- 
sponses through a strong coupling between the fermion parity 
and the spin degree of freedom. The idea behind is schemati- 
cally shown in Fig. 4. The solid lines represent Kitaev chains 
with MBSs indicated by red circles. Let us elaborate a bit 
more on the two kinds of fermions that can be formed by our 
MBSs. Since Majorana fermions are real, normal fermions 
are formed from pairs of them. We can construct two kinds 
of normal (complex) fermions, non-local spin polarized ones 
?/>t = l/2(7n + «7t2) and = 1/2(74,1 - ^2), and lo- 
cal ones ^1 = l/2(7 t i + 274,1) and ip 2 = 1/2(7^2 - ^2). 
The local pseudo spin operators are defined as = 1/2 x 
7i(2)o:< T a/37i(2)/3- It is known that Majorana fermions have 
only Ising-type spin due to their anticommutation relation 
{7i,7j} = 25ij, then, s x = s z = and s y = 1/2 — ip^ip 
from the above definition. We would like to emphasize that 
these local variables are specific to the paired Majorana states. 
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FIG. 3. (a) Density of states with different values of fi. The hori- 
zontal axis is in units of the substrate superconducting gap A. This 
plot shows the topological phase transition at \i — /x* . There are no 
states inside the gap when \i < \i* while there are four-fold degen- 
erate zero energy states when \±> \±* . (b) Probability distribution of 
the zero energy states. The size of system is 160 x 24 for the sub- 
strate and 140 x 1 for the wire. Red solid circles show the probability 
distribution in the wire, which has peaks at both ends. 



There will be various interesting physics as a consequence of 
the synthesis of the local and the non-local nature. Here, as 
an example, we focus on spin-spin correlation of two Ising 
spins formed at the ends. In preparation for that, we first 
review some aspects of the fermion parity. When we have 
two Majorana operators, 7 and 7', there exist two orthog- 
onal states |0) and |1), which satisfy — 277' |0) = |0) and 
— 277' |1) = — 11), respectively. The former has an even 
fermion parity and the latter odd. They can be rewritten by 
normal fermion operators defined as ip = 1/2(7 + 27') where 
^|0) = 0andV^|l) = |1). Note that if the wire is con- 
nected to ground by a capacitor, one can modify the fermion 
number via charging energy U(n) = Q 2 / (2C) — V(t)Q = 
(Q — Qo) 2 /(2C)+const. where Q = ne is the charge, V(t) is 
the gate voltage, and Qo = CV(t). The importance of charg- 
ing energy in the context of topological superconductors was 
first discussed in Ref. [37], where electron teleportation (me- 
diated by the non-locality of the MBSs) has been proposed. 
We now apply this idea to our model. Our Hilbert space is 
4-dimensional and there are two kinds of bases one can con- 
struct, i.e., from ipi, ip 2 , or from and ip^. Namely, one 
can define ^j^ 1 |l)i = |l)i, etc., and |0)i, |l)i, ...are local 
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FIG. 4. Schematic figure of Majorana states at the ends of the wire 
(red circles). As discussed in the main text, the QW can be regarded 
as two copies of Kitaev model (solid lines). There are two ways to 
combine two of them into ordinary fermions: non-local and local 
ones. 



basis states while |0)^, |1)^, . . . are non-local basis states. The 
relation between these two kinds of basis states is given by 
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In our system, the total number of fermions can be controlled 
by the gate voltage V(t). We consider sis 2 correlation in the 
odd and even fermion number space, respectively. The den- 
sity matrices are p odd = \ (|l)i|0) 2 i(0| 2 (l| + |0)i|l) 2 i(l| 2 (0|) 
and/wn = § (|0>i|0> 21 <0| 2 <0| -h | l>i 1 1> 2 i<l 1 2 <1 1 ) . Then we 
obtain 
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where (. . . ) dd/even represents the statistical average with 
the density matrices p Q dd and p even , respectively. These re- 
sults show that by modulating charging energy we can con- 
trol spin-spin correlation, which is independent of the distance 
between the two spins because they are connected through 
the superconducting wire in between. This effect is a con- 
sequence of both, the local and the non-local degrees of free- 
dom of the system. One can also manipulate the two Ising 
spins by modulating the gate voltage as a function of time as 
V(t) = Vo + SV(t). The equations of motion of s\ y and s 2y 
by the driving force SV (t) (iji^^ + £714,724,) yield 



d 
tit 

dH 



(Sly + S 2 y) = 

f2e5V \% , 

[Sly ~ S 2 y) = -4 (— ) {S ly -S 2y ) 



(14) 
(15) 



for a constant 5Vq during some time period T. By tuning T 
to satisfy AeSVoT/h = tt one can reverse the pseudospins 
s y components when they are anti-parallel while they remain 
unchanged when they are parallel. 

Conclusion - In this paper, we have studied theoretically 
the proximity effect of a one-dimensional metallic quantum 
wire without spin-orbit interaction on the substrate of an un- 
conventional superconductor. We have considered three dif- 
ferent cases for the substrate: (i) a chiral superconductor; (ii) 
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a class Dili superconductor with a non-trivial Z 2 number, and 
(iii) a class Dili superconductor with a trivial Z 2 number. In- 
terestingly, we have found degenerate zero energy Majorana 
bound states at both ends of the quantum wire for all three 
cases. We have shown that these degenerate Majorana bound 
states for cases (ii) and (iii) can be nicely used for a combined 
spin/parity qubit. Furthermore, we have demonstrated that the 
resulting non-local spin correlations between the two ends of 
the wire can be controlled by the gate voltage potential acting 
on the wire. One of the candidate systems for the case (ii) 
substrate are thin films of Sr2Ru04 [38] or bi-layer Rashba 
system [30], where the class Dili superconducting state might 
be realized, and for the case (iii) super lattice structure of 
CeCoIn 5 [39]. 
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